Full Numerical Estimation of Neutrino Mixing Parameters from Solar Neutrino Data 
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We have numerically computed survival probabilities of solar neutrinos interacting with matter 
via MSW mechanism in the full three generation formalism with Am^ = 2.2 x 10~ 3 eV 2 . For 
#13 < 30°, we confirm the two regions found in the two generation case, (1) with the most likelyhood, 
small sin 2 (20i 2 ) ~ 0.006, (2) with less likelyhood, large sin 2 (26»i 2 ) > 0.5. For 0i 3 > 30°, we find 
• an additional region, (3) with even less likelyhood, 10" 5 < Am 2 2 < 10 -4 eV 2 and/or 10~ 4 < 

| sin 2 (26>i 2 ) < 0.5. Assuming that the case (1) is valid, we predict that #13 ~ 15°. 

PACS numbers: 14.60.Pq; 26.65.+t; 96.60. Jw 

Recent Super-Kamiokande experiment jjj shows an evidence for oscillations of the atmospheric neutrinos. The data 
are in good agreement with two-flavor <-> v T oscillations. These results did establish that neutrinos oscillate and 
possess non-zero masses. Thus full three-flavor neutrino oscillation analysis is needed for more accurate determination 
of neutrino masses and mixing anlges. There are many investigations of neutrino oscillations with three neutrino 
generation schemes of solar neutrinos, accelerator and reactor neutrinos as well as atmospheric neutrinos |^ |J. 
The solar neutrino puzzle is the discrepancy between the standard solar model predictions jp](| and the results of 
qq | the solar neutrino experiments. The standard solar model of Bahcall and Pinsonneault predicts that the neutrino 
£N) • event rates are 9.3^' 2 SNU for the chlorine experiment, 137^y SNU for the gallium experiment, and (6.62^i fi) x 10 6 
/cm 2 /sec for the Super-Kamiokande experiment. The observed neutrino event rates are 2.56 ± 0.16 ± 0.15 SNU for 
Homestake Chlorine experiment [§, 77.5 ± 6.2±|;f SNU for GALLEX Gallium experiment g, 67.21^1^ SNU for 
SAGE Gallium experiment H), and (2.44 ± 0.05too?) x 106 /cm 2 /scc for the Kamiokande experiment |To| . 
■ Several mechanisms were proposed for the depiction of the solar neutrinos. Among them the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect |fj"l[ , [l2| has been most popular because it induces large flavor transitions in spite of small 
vacuum mixing angle in an energy dependent way. Two-flavor MSW solutions |L3|Jl4j| to the solar neutrino puzzle was 
I \ fully analysed with high accuracy, including both the day/night effect and spectral information. 

In the two-neutrino scheme with the MSW effect it was found that there are two regions in the sin 2 29 — Am 2 
parameter space that yield two different solutions of the solar neutrino problem: the small mixing angle solution and 
the large mixing angle solution. The chi-square analysis indicates that the small mixing angle solution is favored p3[. 
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Authors in ref. |15| showed that the high- Am 2 part of the large mixing angle MSW solution of the solar neutrino 



problem is disfavored by the Super-Kamiokande atmospheric neutrino data assuming the scheme of neutrino mixing 



, — i 

x. 

indicated by the result of the reactor neutrino oscillation experiment CHOOZ |16|]. This implies that the mixing 
matrix element U e 3 (or #13 in our parametrization which is the same as that of the Particle Data Group ) is small 
and the oscillations of solar and atmospheric neutrinos are decoupled under the assumption of neutrino mass hierarchy 
Am 2 2 <C Amf.3, (Ato 2 2 = to 2 . — to 2 , Am| 3 = to 2 — m|). 

Since the solar neutrino oscillates into two other flavors, a full analysis may result in some new findings. Especially 
the angle #13 may play a non-trivial role. There are many papers in which the MSW effect is considered in three- 
generation-neutrino scheme 18-2^]. However, they used the Landau-Zenor approximate formula for treating flavor 



conversion in the non-adiabatic resonant region of the solar interior with the assumption that all neutrinos are created 



at the center of the Sun. Two groups [|26 27 implemented numerical methods with two neutrino flavors. 

This method of using the Landau-Zenor formula is unreliable when the resonant region and the neutrino production 
region overlap. In order to demonstrate the situation, we first present a sample of our numerical analysis which shows 
that for low energy neutrinos the survival probabilities indeed depend very much on their creation sites, r. The survival 
probability of a solar neutrino with energies, 0.3 < E v < 1.0 MeV, varies rapidly in the range, < r/i? Q < 0.2. Thus 
r-dependence of the survival probability should be taken into account for the pp and Be neutrinos though not as much 
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for the boron neutrinos. That is, one should consider the radial distribution of the solar neutrinos in analysing the 
gallium experiment data [p|j9|j2^] and the Borexino data |2fi| . 

Recently we [ j30| presented a numerical algorithm for computing the survival probability of an electron neutrino 
in its flight through the solar core experiencing the MSW effect. We adopted a hybrid method to integrate the 
neutrino evolution equation in the full three-generation formalism. In order to reduce the number of integrations we 
used the importance sampling method for sampling the neutrino creation energy and position. In order to reduce 
the amount of computation in each integration we need to find the optimum radii at which numerical integration is 
started and stopped. Outside these radii we would use the adiabatic conversion formula. We checked the location of 
the non-adiabatic resonant regions, where direct numerical integration is performed, by checking if the relative ratios 
of the off-diagonal and diagonal elements of the Hamiltonian are below a certain limit, say j c . Thus the frequency 
and amount of numerical integrations are greatly reduced. Even with these tactics it is extremely time consuming to 
solve the neutrino evolution equations for the three-generation case in a large region of the parameter space and it 
requires deployment of a parallel supercomputer. We developed a parallel algorithm for a message passing parallel 
computer. 

In this paper, we present results of our lengthy numerical computation of the solar neutrino survival probabilities. 
We will show that the numerical MSW solutions to the solar neutrino problem within the three-generation formalism 
lead to a new constraint to the neutrino masses and mixing angles. In the three neutrino formalism, we have three 
vacuum mixing angles (9i2, O23, #13) and one phase to describe mixing between mass and flavor eigenstates. The 
angle #23 is irrelevant in the solar neutrino problem and the CP violating phase may be ignored. Thus adjustable 
parameters are two mixing angles and two mass squared differences. In conformity with the Super-Kamiokande data 
on the atmospheric neutrinos, we fixed two parameters, Am^ = 2.2 x 10 -3 eV 2 and 823 = 43.5° and searched a 
wide range of the parameter space spanned by Am 2 2 , 6*12, and #13. We used a 24x21 logarithmic grid for Am 2 2 and 
sin 2 (26 l i2). We focused on the interesting region, 2 x 1CT 6 < Am 2 2 < 4 x 1CT 4 eV 2 and 1CT 4 < sin 2 (26»i 2 ) < 1. 

In our production run, we refined the importance sampling method for sampling neutrino creation energies and 
radii. For the boron and CNO neutrinos, we chose 64 energy sample points ranging from 0.23 MeV to 18.65 MeV and 
60 radius sample points ranging from to 0.353i?Q. The three energy values, 0.384, 0.862, 1.442 MeV, are included to 
evaluate the Be and pep neutrinos more accurately. For the most numerous and model independent pp neutrinos, to 
be counted by the gallium detectors, we chose 60 energy sample points ranging from 0.23 MeV to 0.4293 MeV and 36 
radius sample points ranging from to 0.4i? Q . In order to determine the radial limits of the non-adiabatic region, we 
used a very generous critical value 7 C = 0.005, which guaranteed the accuracy of 0.001 for the survival probabilities 
in the worst case. 

We did not include the Earth effect in the calculation. The sensitive parameter region to the Earth density is 
around Am 2 2 ~ 10 -6 eV 2 . In the two-generation scheme the Earth effect does not cause much change of the 
combined parameter regions for the solar neutrino problems. The data of energy spectrum and day/night variation 
from the Kamiokande and Super-Kamiokande may give more restriction on the parameter space of neutrino mass and 
mixing angle. The absence of any day/night variation in the data is consistent with neglecting the MSW resonance 
effect in the Earth for the electron neutrino. The excluded region by the data of neutrino energy spectrum is around 
Am 2 2 ~ 10~ 4 eV 2 which is outside the combined parameter regions in the two-generation case. 

To compute the event rates, we used the 98 SSM flux data taken from (5^]. We present our results for several 
values of # 13 in the familiar iso-SNU/FLUX contour format in Figs. 1. In these figures only the mean SNU(FLUX) 
values are plotted. From these figures we notice three important facts: (1) When <9 13 < 22.5°, the small mixing angle 
region where sin 2 (26' 12 ) <C 1, is the only region where all three curves can come close. In the large mixing angle 
regions two curves can cross at a time but the crossing points are rather far apart. (2) Beyond #13 > 22.5°, the gap 
which the Ga curve passes between the CI curve and the Water curve becomes wide and even the small mixing angle 
solution becomes less likely. (3) When #13 > 35°, there is no place in the considered region of the Am 2 2 — sin 2 (2#i 2 ) 
plane where any pair of curves can meet. Thus the likelyhood to have a solution with a large (1-3) mixing angle is 
very small. 

As we can see in the figures the chlorine data are least sensitive to 613. The Water data are a little more sensitive 
to 6>i3 but the shape of the curve does not change. In the small Am 2 2 region with small values of #13 the Water 
curves overlap the CI curves as both experiments are sensitive mostly to the boron neutrino. Thus it is impossible 
to determine the mixing angles decisively from these two experiments. However, the gallium detector is sensitive to 
the pp neutrinos also and its iso-SNU curve crosses the other two curves at large angles and thus most useful for a 
precision determination of #i 2 . The figures show that the gallium curves are also very sensitive to a choice of #13 and 
they can be used to impose a non-trivial constraint on #13. 

Experimental data contain two kinds of errors, statistical and systematic. We take weighted averages and estimate 
combined errors as reviewed in jlTj. For the CI experiment, we use Rci = 2.56 ± 0.23 SNU and for the Water 
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experiment, we used i?water = (2.44±0.10) x 10 6 /cm 2 /sec. For the two gallium experiments we take weighted average 
to get i?Ga = 73.1 ± 6.6 SNU. The flux data from the SSM model contain some theoretical uncertainties. Thus we 
need to represent the la deviation values of SNU(FLUX). Figs. 2 take into account both theoretical uncertainties 
and experimental errors. The upper curves are from the lower limit of the SSM flux data and the upper values of 
detection rates of each experiment and vice versa. 

For small values of 0i3, our result is consistent with the two generation case Jl^,[l4| as expected. Barbieri et. al. 
p5[ performed extensive analysis fo solar neutrino experiments with all three neutrinos. They used the Landau-Zenor 
approximation formula given in Jl^JTsj ] - The experimental data are slightly updated since their work. They used, 
Ra = 2.54 ± 0.20 SNU for CI, i? Ga = 75 ± 7 SNU for Ga, and fl W ater = (2.51 ± 0.16) x 10 6 /cm 2 /sec for Water. So 
direct comparison of our results with their results is not possible. However, we made iso-SNU(FLUX) plots with these 
data also. For #13 = 1.5°, 15°, our best estimates of Ato 2 2 in units of 10~ 6 eV 2 and sin 2 (20i 2 ) in units of 10 -3 are 
(5.05, 5.89) and (6.50, 3.99) respectively. However, for 6*13 = 30° our estimates are (7.44, 2.20), whereas their plots 
indicate (~ 8, ~ 6). It seems that their estimates of Am 2 2 are slightly larger than ours but still within the la limit. 
For small values of 813 we agree with theirs within the la limit whereas for large values our estimate of sin 2 (20i 2 ) is 
much smaller. In both cases, for large #13, there is a la region where Am 2 2 ~ 10~ 4 eV 2 . 

Due to the theoretical uncertainties and experimental errors, all angular values of #13 can accomodate all three 
experiments in broad ranges of Am 2 2 and Qyi- For small #13 the 'small mixing angle' solution gives the most likely 
solution, with a small patch of 'large mixing angle' region hanging in the upper-right corner. However, for large #13, 
the large boomerang shaped region in Figs. 2 allows broad range of values for Am 2 2 and sin 2 (26*i 2 ). Large values of 
Am 2 2 w 10 -4 eV 2 are possible and 2 • 10~ 4 < sin 2 (26>i 2 ) < 0.2 for 6>i 3 = 40°. But its likelyhood seems to be minimal 
as Figs. 1 indicate. However, the earth effect can broaden the large 812 region, which we did not consider in this 
work. 

We repeated the computation with Am 23 = 5.0 x 10 -4 eV 2 , which is the la lower limit of the Super-Kamiokande 
estimates. For small values of 6*13 < 22.5° the contours of mean SNU values (Fig. 1) look almost the same as in 
the case Am 23 = 2.2 x 10~ 3 eV 2 . However, in the latter case (Am| 3 = 5.0 x 10 -4 eV 2 ), for large values of #13 the 
Water contours are shifted slightly to the left and the Ga contour are pushed further to the upper-left direction, while 
the CI contours are intact. Nevertheless we do not see noticablc differences between the two cases. However, with 
Am 23 = 2.0 x 10~ 4 eV 2 , we noticed non-trivial deviations. In this case (2-3) transitions can take place in addition 
to the (1-2) transition. Our program has been improved to compute the survival probabilities even for the most 
general situation when all three transitions occur. We repeated the computation with the same set of 6*13. The most 
significant difference is that the la region broadens for large 613. 

We give the values of the small (1-2) mixing angle solutions within la range in Table 1. It is interesting to notice 
that Xmin depends on #13. Its minimum value occurs when the CI and the Water curves overlap perfectly. It can 
be used to predict the most likely value of #13 from solar neutrino data alone. In this case the minimum occurs at 
013 = 13.21°. 



Table 1. Best estimates of small (1-2) mixing parameters for several values of $13 obtained by minimizing 
X 2 . Here, we used the definition x 2 = Ylid^-i ~ R° xp ) / AR°* P ) 2 with i = CI, Ga, Water, which is caused 
purely by the experiments. The la bounds are also indicated. 
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Xmin 



6.29: 



2.81] 
2.12 H 



0.04526 
0.05902 
0.07068 
0.00237 
0.03779 
0.76387 
3.36182 



In summary we have computed the solar neutrino survival probabilities P(v e — > 1/3; E v ,r] Am?-,0y) (1) within the 
full three generation framework (2) using accurate numerical methods and (3) included the effects due to neutrino 
creation positions to full extent. We have found that there are three regions in the parameter space where all three 
solar neutrino experiments can be accomodated within the la limit: for small (1-3) mixing angles 6*13 < 30°, (1) 
with the most likelyhood, small A?ti 2 2 ~ 5 x 10 -6 eV 2 and small sin 2 (20i 2 ) ~ 0.006, (2) with less likelyhood, large 
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Amj 2 ~ 10 4 eV 2 and large sin 2 (26*12) > 0.5; and for large (1-3) mixing angles #13 > 30°, (3) with even less likelyhood, 
10- 5 < Am 2 2 < 10~ 4 eV 2 and/or 1CT 4 < sin 2 (26» 12 ) < 0.5. 

From Figs. 1 and Figs. 2 and Table 1, we conclude that the small sin 2 (26*12) solution with 6*13 smaller than 22.5° 
will cope with future experiments more likely than the large sin 2 (26*12) solutions. We estimate the most likely value 
of 6>i3 by minimizing Xmm w ^ n respect to 6 , i 3 , which comes out to be 6*13 = 13.21° from the Table 1. Our conclusion 
is consistent with the CHOOZ estimate of 6*13 < 15°. If the errors in the CI and the Water experiments are reduced 
further and the value of 6*13 measured in CHOOZ and other experiments agree with our prediction, it can serve as an 
interesting test of the Standard Solar Model and the MSW mechanism. 

We would like to thank Dr. J.D. Kim for discussions and Ms. K.S. Lee for making plots and analysing data for 
the table. We also thank Dr. S.B. Kim for clarifying the Super-Kamiokande data. This work was funded by the 
POSTECH research fund. 
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Fig. 1. iso-SNU(FLUX) plots for solar neutrinos with the CI, Ga, and Water detectors. 

Contour values are 2.56 SNU (CI), 73.1 SNU (Ga) and 2.44x 10 6 /cm 2 /sec (Water) and they are represented 
by solid, dash-dot, and dotted curves respectively. In each plot the vertical axis refers to Am 2 2 in eV 2 
and the horizontal axis to sin 2 (2#i 2 ). We used Am| 3 = 2.2 x 10~ 3 eV 2 . 
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